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We review recent advances in the theory of strong-interaction effects and final-state interactions in hadronic 
weak decays of heavy mesons. In the heavy-quark limit, the amplitudes for most nonleptonic, two-body B 
decays can be calculated from first principles and expressed in terms of semileptonic form factors and light- 
cone distribution amplitudes. We summarize the main features of this novel QCD factorization and discuss its 
phenomenological applications to B — > £)'*' L decays (with L a light meson), and to the rare charmless decays 
B —* nK and B — > tttv. 



1. INTRODUCTION 

The theoretical description of hadronic weak 
decays is difficult due to nonperturbative strong- 
interaction dynamics. This affects the interpreta- 
tion of data collected at the B factories, and in 
many cases limits our ability to uncover the origin 
of CP violation and search for New Physics. The 
complexity of the problem is illustrated in Fig. [l| 

It is well known how to control the effects 
of hard gluons with virtuality between the elec- 
troweak scale Mw an d the scale tub characteris- 
tic to the decays of interest. They can be dealt 
with by constructing a low-energy effective weak 
Hamiltonian 

H ^ = ^ E ^ C i( M w/») Oi(jM) + h.C, (1) 
* i 

where Ap KM are products of CKM matrix ele- 
ments, Ci(Mw I V) are calculable short-distance 
coefficients, and Oi(fi) are local operators renor- 
malized at a scale /i = 0(ms). The challenge is 
to calculate the hadronic matrix elements of these 
operators with controlled theoretical uncertain- 
ties, using a systematic approximation scheme. 

Previous field-theoretic attempts to evaluate 
these matrix elements have employed dynami- 
cal schemes such as lattice field theory, QCD 
sum rules, or the hard-scattering approach. The 
first two have difficulties in accounting for final- 
state rescattering, which however is important 
for predicting direct CP asymmetries. The hard- 
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Figure 1. QCD effects in a hadronic decay. 

scattering approach misses the leading soft con- 
tribution to the B — > meson transition form fac- 
tors and thus falls short of reproducing the cor- 
rect magnitude of the decay amplitudes. In view 
of these difficulties, most previous analyses of 
hadronic decays have employed phenomenologi- 
cal models such as "naive" or "generalized fac- 
torization" , in which the complicated matrix el- 
ements of four-quark operators in the effective 
weak Hamiltonian are replaced, in an ad hoc way, 
by products of current matrix elements. Correc- 
tions to this approximation are accounted for by 
introducing a set of phenomenological parameters 
a,i. A different strategy is to classify nonleptonic 
decay amplitudes according to flavor topologies 
("trees" and "penguins"), which can be decom- 
posed into SU(3) or isospin amplitudes. This 
leads to relations between decay amplitudes in 
the flavor-symmetry limit. No attempt is made 
to compute these amplitudes from first principles. 
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2. QCD FACTORIZATION FORMULA 

Here we summarize recent progress in the the- 
oretical understanding of nonleptonic decay am- 
plitudes in the heavy-quark limit PL0]. The un- 
derlying idea is to exploit the presence of a large 
scale, i.e., the fact that » Aqcd- In order 
to disentangle the physics associated with these 
two scales, we factorize and compute hard con- 
tributions to the decay amplitudes arising from 
gluons with virtuality of order mb, and parame- 
terize soft and collinear contributions. Consider- 
ing the cartoon in Fig. [l], we denote by Mi the 
meson that absorbs the spectator quark of the 
B meson, and by M2 the meson at the upper 
vertex, to which we refer as the "emission parti- 
cle" . We find that at leading power in Aqcd /mb 
all long-distance contributions to the decay am- 
plitudes can be factorized into semileptonic form 
factors and meson light-cone distribution ampli- 
tudes, which are much simpler quantities than the 
nonleptonic amplitudes themselves. A graphical 
representation of the resulting "factorization for- 
mula" is shown in Fig. |^. The physical picture un- 
derlying factorization is color transparency HQ. 
If the emission particle is a light meson, its con- 
stituents carry large energy of order mb and are 
nearly collinear. Soft gluons coupling to this sys- 
tem see only its net zero color charge and hence 
decouple. Interactions with the color dipole of the 
small gg-pair are power suppressed in the heavy- 
quark limit. 

For B decays into final states containing a 
heavy charm meson and a light meson, the fac- 
torization formula takes the form 

(D^ + L-\O l \B d )=J2Ff- Di ' ) f L 

3 

1 

x JduT£(u)$ L (u) + o(^y (2) 


where Oi is an operator in the effective weak 
Hamiltonian ([!]), fB^d { ' are transition form 
factors (evaluated at q 2 = m\ ~ 0), /l and <I>l 
are the decay constant and leading-twist light- 
cone distribution amplitude of the light meson, 
and T-j are process-dependent hard-scattering 
kernels. For decays into final states containing 




Figure 2. QCD factorization in the heavy-quark 
limit. The second term is power suppressed for 
B — > Dn, but must be kept for decays with two 
light mesons in the final state, such as B — > irK. 
Contributions not shown (such as weak annihila- 
tion graphs) are power suppressed. 

two light mesons there is a second type of con- 
tribution to the factorization formula, which in- 
volves a hard interaction with the spectator quark 
in the B meson. This is shown by the sec- 
ond graph in Fig. [|. Below we focus first on 
B — > D^L decays (with L a light meson), where 
this second term is power suppressed and can 
be neglected. Decays into two light final-state 
mesons are more complicated jl],|| and will be 
discussed briefly in Sect. ||. 

The factorization formula for nonleptonic de- 
cays provides a model-independent basis for the 
analysis of these processes in an expansion in 
powers and logarithms of Aqcd /mb- At leading 
power, but to all orders in a s , the decay ampli- 
tudes assume the factorized form shown in (^). 
Having such a formalism based on power count- 
ing in Aqcd /mb is of great importance to the 
theoretical description of hadronic weak decays, 
since it provides a well-defined limit of QCD in 
which these processes admit a rigorous, theoret- 
ical description. (For instance, the possibility 
to compute systematically 0{a s ) corrections to 
"naive factorization" , which emerges as the lead- 
ing term in the heavy-quark limit, solves the old 
problem of renormalization-scale and scheme de- 
pendences of nonleptonic amplitudes.) The use- 
fulness of this new scheme may be compared with 
the usefulness of the heavy-quark effective theory 
for the analysis of exclusive semileptonic decays 
of heavy mesons, or of the heavy-quark expan- 
sion for the analysis of inclusive decay rates. In 
all three cases, it is the fact that hadronic uncer- 
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tainties can be eliminated up to power corrections 
in Aqcd / mb that has advanced our ability to con- 
trol theoretical errors. 

It must be stressed, however, that we are just 
beginning to explore the theory of nonleptonic B 
decays. Some important conceptual problems re- 
main to be better understood. In the next few 
years it will be important to further develop the 
approach. This should include an all-orders proof 
of factorization at leading power, the develop- 
ment of a formalism for dealing with power cor- 
rections to factorization, understanding the light- 
cone structure of heavy mesons, and understand- 
ing the relevance (or irrelevance) of Sudakov form 
factors. Also, we must gauge the accuracy of the 
approach by learning about the magnitude of cor- 
rections to the heavy-quark limit from extensive 
comparisons of theoretical predictions with data. 

As experience with previous heavy-quark ex- 
pansions has shown, this is going to be a long 
route. Yet, already we have obtained important 
insights. Before turning to specific applications, 
let me mention three points here: 

1. Corrections to "naive factorization" (usu- 
ally called "nonfactorizable effects") are process 
dependent, in contrast with a basic assumption 
underlying models of "generalized factorization" . 

2. The physics of nonleptonic decays is both 
rich and complicated. There may, in general, be 
an interplay of several small parameters (Wil- 
son coefficients, CKM factors, 1/N C , etc.) in 
addition to the small parameter AQcrj/mf, rel- 
evant to QCD factorization. Also, several not- 
so-well- known input parameters (e.g., heavy- to- 
light form factors and light-cone distribution am- 
plitudes) introduce numerical uncertainties. 

3. Strong rescattering phases arising from final- 
state interactions are suppressed in the heavy- 
quark limit. More precisely, the imaginary parts 
of nonleptonic decay amplitudes are suppressed 
by at least one power of a s {mb) or Aqcd At 
leading power, the phases are calculable from the 
imaginary parts of the hard-scattering kernels in 
the factorization formula. A generic consequence 
is that direct CP asymmetries in hadronic B de- 
cays are suppressed in the heavy-quark limit, ex- 
cept for cases where for some reason the real part 
of the decay amplitude is also suppressed. 



3. APPLICATIONS TO B d -> DW+i" 
DECAYS 

Our result for the nonleptonic Bd — > D^ + L^ 
decay amplitudes (with L a light meson) can be 
compactly expressed in terms of the matrix ele- 
ments of a "transition operator" 

Gp 



T 



v/2 



v: d v cb 



ai {DL)Q v - ai {D*L)Q A , (3) 



where the hadronic matrix elements of the op- 
erators Qv = ey^b ® <?7/j(1 — 75) u an d Qa = 
C7^7 5 & (g> (£7^(1 — 75 )u are understood to be evalu- 
ated in factorized form. Eq. ([}]) defines the quan- 
tities ai(D^L), which include the leading "non- 
factorizable" corrections, in a renormalization- 
group invariant way. To leading power in 
Aqcd/w& these quantities should not be inter- 
preted as phenomenological parameters (as is 
usually done), because they are dominated by 
hard gluon exchange and thus calculable in QCD. 
At next-to-leading order in a s we obtain B 



/ n {.) rN f=t I \ I ^( m b) 

ai(D ( >L) = Ci(m b ) -I — — 



1 



47T 



X 

J duF(u,±z) $ L (u) 



, (4) 



where z — nic/mi,, Ciirrib) are the so-called 
"renormalization-scheme independent" Wilson 
coefficients, and the upper (lower) sign in the sec- 
ond argument of the function F(u, ±z) refers to 
a D (D*) meson in the final state. The exact ex- 
pression for this function is known but not rele- 
vant here. Note that the coefficients a\(DL) and 
a\(D*L) are nonuniversal, i.e., they are explic- 
itly dependent on the nature of the final-state 
mesons. Politzer and Wise have computed the 
"nonfactorizable" vertex corrections to the decay 
rate ratio of B d -> D+n- and B d -> D*+ix- §. 
This requires the symmetric part (with respect to 
u <-> 1 — u) of the difference F(u, z) — F(u, —z). 
We agree with their result. 

The expressions for the decay amplitudes ob- 
tained by evaluating the hadronic matrix ele- 
ments of the transition operator T involve prod- 
ucts of CKM matrix elements, light-meson de- 
cay constants, B — > transition form factors, 
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and the QCD parameters a\{D^*'L). A numer- 
ical analysis shows that |oi| = 1.055 ± 0.025 for 
the decays considered below Below we will 
use this as our central value. 

3.1. Tests of factorization 

A particularly clean test of our predictions 
is obtained by relating the Bd — > D* + L~ de- 
cay rates to the differential semileptonic Bd — > 
D* + l~v decay rate evaluated at q 2 — m\. In 
this way the parameters |oi| can be measured di- 
rectly, since 

T(B d -> D*+L-) 
dT{B d ^D*+l-D)/dq>\ q2=ml 

= ^ 2 \V ud \ 2 f L \a 1 (D*L)\ 2 . (5) 

With our result for a\ this relation becomes a 
prediction based on first principles of QCD. This 
is to be contrasted with the usual interpretation 
of this formula, where a\ plays the role of a phe- 
nomenological parameter that is fitted from data. 

Using results reported by the CLEO Collabo- 
ration |7j], we find 

\ ai {D*Tt)\ = 1.08 ±0.07, 
\ ai {D*p)\ = 1.09 ±0.10, 

\a x {D*ai)\ = 1.08 ±0.11, (6) 

in good agreement with our prediction. It is re- 
assuring that the data show no evidence for large 
power corrections to our results. However, a fur- 
ther improvement in the experimental accuracy 
would be desirable in order to become sensitive 
to process-dependent, nonfactorizable effects. 

3.2. Predictions for class-I amplitudes 

We now consider a larger set of so-called class-I 
decays of the form Bd — > Z?'*' + L~, all of which 
are governed by the transition operator (|J). In 
Tab. |] we compare the QCD factorization pre- 
dictions with experimental data. As previously 
we work in the heavy-quark limit, i.e., our pre- 
dictions are model independent up to corrections 
suppressed by at least one power of h.QCT> /rrib- 
There is good agreement between our predictions 
and the data within experimental errors, which 
however are still large. It would be desirable to re- 
duce these errors to the percent level. (Note that 



Table 1 

Model-independent predictions for the branching 
ratios (in units of 10" 3 x (| ffll |/1.05) 2 ) of B d -> 
£)(*)+ L~ decays in the heavy-quark limit. 



Decay mode 


Theory (HQL) 


PDG § 


B d -> D+-R- 
B d -► D+K- 
B d D+p- 
B d -» D+K*- 
B d -> D+a^ 


3.27 x [F+/0.6] 2 
0.25 x [F+/0.6] 2 
7.64 x [F+/0.6] 2 
0.39 x [F+/0.6] 2 
7.76 x [.F+/0.6] 2 


3.0 ±0.4 
7.9 ± 1.4 
6.0 ±3.3 


B d D*+n- 
B d -> D*+K~ 
B d D*+p- 
B d -► D*+K*- 
B d D*+ai 


3.05 x [Aq/0.6] 2 
0.22 x Uo/0.6] 2 
7.59 x [Ao/0.6] 2 
0.40 x Uo/0.6] 2 
8.53 x Uo/0.6] 2 


2.8 ±0.2 
6.7 ±3.3 
13.0 ±2.7 



we have not attempted to adjust the semileptonic 
form factors F+^ D and A®^ D entering our re- 
sults so as to obtain a best fit to the data.) 

The observation that the experimental data on 
class-I decays into heavy-light final states show 
good agreement with our predictions may be 
taken as (circumstantial) evidence that in these 
decays there are no unexpectedly large power cor- 
rections. In our recent work [|| we have addressed 
the important question of power corrections theo- 
retically by providing estimates for two sources of 
power-suppressed effects: weak annihilation and 
spectator interactions. A complete account of 
power corrections to the heavy-quark limit can 
at present not be performed in a systematic way, 
since these effects are no longer dominated by 
hard gluon exchange. However, we believe that 
our estimates are nevertheless instructive. 

We parameterize the annihilation contribution 
to the Bd — > D + n~ decay amplitude in terms 
of an amplitude A such that A(Bd — * D + u~ ) = 
T + A, where T is the "tree topology" , which 
contains the dominant factorizable contribution. 
We find that A/T - 0.04. We have also ob- 
tained an estimate of nonfactorizable spectator 
interactions, which are part of T, finding that 
Tgpec/Xiead ~ —0.03. In both cases, the results 
exhibit the expected linear power suppression in 
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the heavy-quark limit. We conclude that the typ- 
ical size of power corrections in class-I decays into 
heavy-light final states is at the level of 10% or 
less, and thus our predictions for the values and 
the near universality of the parameters a\ gov- 
erning these decay modes appear robust. 

4. QCD FACTORIZATION IN CHARM- 
LESS HADRONIC B DECAYS 

The observation of rare charmless B decays 
into nK and tttt final states has resulted in a large 
amount of theoretical and phenomenological work 
that attempts to interpret these observations in 
terms of the factorization approximation, or in 
terms of general parameterizations of the decay 
amplitudes. A detailed understanding of these 
amplitudes would help us to pin down the value of 
the CKM angle 7 = axg(V* b ) using only data on 
CP-averaged branching fractions. Here we briefly 
summarize the most important consequences of 
the QCD factorization approach for the irK and 
7T7T final states B . 

To leading order in an expansion in powers of 
A-QCD/^fc, the B — > ttK matrix elements obey 
the factorization formula shown in Fig. [2J 

(irK\0 i \B)=F*^ r f K Tj Cti *9 K 

+ /B/^T/ 7 *$ fl *$ K *$ 7r , (7) 

where the ^-products imply an integration over 
the light-cone momentum fractions of the con- 
stituent quarks inside the mesons. Our results 
are based on hard-scattering kernels including all 
corrections of order a s . Compared to our previ- 
ous discussion of B — > nn decays [Q, the present 
analysis incorporates three new ingredients: the 
matrix elements of electroweak penguin operators 
(for ttK modes) , hard-scattering kernels for asym- 
metric light-cone distributions, and the complete 
set of "chirally enhanced" corrections. The 

second and third items have not been considered 
in other generalizations of Ref. jy to the ttK final 
states |],[l(J. The third one, in particular, is es- 
sential for estimating some of the theoretical un- 
certainties of the approach. For completeness, we 
note that the predictions from QCD factorization 



differ in essential aspects from those obtained in 
the conventional hard-scattering approach Jljj . 

Following Ref. §J, we have obtained the co- 
efficients ai(irK) (with i = 1, ...,10) of the ef- 
fective, factorized "transition operator" defined 
analogously to the case of B — > nn decays, but 
augmented by coefficients related to electroweak 
penguin contributions. Chirally enhanced cor- 
rections arise from twist-3 two-particle light-cone 
distribution amplitudes, whose normalization in- 
volves the quark condensate. The relevant pa- 
rameter, 2jjL- K /m,b = —4(qq)/{f„mb), is formally 
of order Aqcd /mi, but large numerically. The co- 
efficients ag and ag are multiplied by this param- 
eter. There are also additional chirally enhanced 
corrections to the spectator-interaction term in 
, which turn out to be the more important 
effect. These latter corrections involve a loga- 
rithmically divergent integral, which violates fac- 
torization. For instance, for matrix elements of 
V — A operators the hard spectator interaction is 
proportional to (u = 1 — u, v = 1 — v) 
1 

/du dv , , . / . . 2/1^ u\ 
— — **(») + — - 8 

u v \ m b uj 



when the spectator quark goes to the pion. (Here 
we used that the twist-3 distribution amplitudes 
can be taken to be the asymptotic ones when 
one neglects twist-3 corrections without the chiral 
enhancement.) The divergence of the ^-integral 
in the second term as v — ► implies that it is 
dominated by soft gluon exchange between the 
spectator quark and the quarks that form the 
kaon. We therefore treat the divergent integral 
X = J^(dv/v) as an unknown parameter, which 
may in principle be complex owing to soft rescat- 
tering in higher orders. In our numerical analy- 
sis we set X — ln(ms/0.35 GeV) + r, where r is 
chosen randomly inside a circle in the complex 
plane of radius 3 ("realistic") or 6 ("conserva- 
tive"). Our results also depend on the B-meson 
parameter [jj] As, which we vary between 0.2 and 
0.5 GeV. Finally, there is in some cases a nonneg- 
ligible dependence of the coefficients ^(ttK) on 
the renormalization scale, which we vary between 
nib/ 2 and 2mb- 

We take \V ub /V cb \ = 0.085 and m s (2GeV) = 
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Figure 3. Ratios of CP-averaged B — > irK and irir decay rates. The scattered points cover a realistic 
(dark) and conservative (light) variation of input parameters. The dashed curve is the result obtained 
using "naive factorization" . 



HOMeV as fixed input to our analysis, noting 
that ultimately these Standard Model parame- 
ters, along with the CP-violating phase 7, might 
be extracted from a simultaneous fit to the B — ► 
ttK and B — > ttt: decay rates. We now summarize 
our main results. 

4.1. Results on SU(3) breaking 

Bounds on 7 derived from ratios of CP- 
averaged B — > ttK decay rates Jl^,|l^|, as well 
as the determination of 7 using the method of 
Ref. |0, rely on an estimate of SU(3) flavor- 
symmetry violations. We find that "nonfactoriz- 
able" SU(3)-breaking effects (i.e., effects not ac- 
counted for by the different decay constants and 
form factors of pions and kaons in the conven- 
tional factorization approximation) do not exceed 
the level of a few percent. 

4.2. Ratios of CP-averaged rates 

The approach discussed here allows us to cal- 
culate the amplitudes for B — ► irK and B — > irir 
decays in terms of form factors and light-cone 



distribution amplitudes. We focus on decays 
whose branching ratios have already been mea- 
sured. Since the relevant form factor .F+ _>7r (0) is 
not well known, we consider only ratios of CP- 
averaged branching ratios. We display these as 
functions of the weak phase 7 in Fig. Ej. For com- 
parison, we show the data on the various ratios 
obtained using results reported by the CLEO Col- 
laboration E3j. We also indicate the very recent 
results reported by the BaBar and Belle Collabo- 
rations at the ICHEP2000 Conference @. From 
the plots it is apparent that the corrections with 
respect to the conventional factorization approx- 
imation are significant (and important to reduce 
the renormalization-scale dependence). Despite 
this fact, we find that the qualitative pattern 
emerging for the set of B — > ttK and tttt decay 
modes is similar to that in conventional factor- 
ization. In particular, the penguin-tree interfer- 
ence is constructive (destructive) in B — > tt + tt~ 
(B — > tt~K + ) decays if 7 < 90°. Taking the cur- 
rently favored range 7 = (60 ± 20)°, we find the 
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following robust predictions: 
Br(7r+7r~) 

4^4=0.5-1.9 [0.25 ±0.10] 

^A__i = 0.9-1.3 [1.27 ±0.47] 

— l^fSnl = 0-6-1-0 [1-00 ±0.30] 

Br( f f = 0-9-1-4 [0.59 ±0.27] 

2Br(7r°if°) 1 J 



sin(2P)=0.75 



(9) 



The first ratio is clearly in disagreement with cur- 
rent CLEO data |lj] shown in square brackets. 
However, there is good agreement with the re- 
cent results reported by BaBar (0.74 ± 0.29) and 
Belle (0.36 ±0.26) @. 

The near equality of the second and fourth ra- 
tios in (|9|) is a consequence of isospin symmetry 
13). We find Br(B -> tt°K°) = (4.5 ± 2.5) x 
6 x [Ff - > ' r (0)/0.3] 2 almost independently of 
This is three time smaller than the central 



10 
7- 



value reported by CLEO, (14.6 



+5.9+2.4 
-5.1-3.3 



) X 10" 



and four times smaller than the central value re- 



ported by BELLE, (21.0 



+9.3+2.5\ 
-7.8-2.3J 



X 10" 



It will 



be interesting to follow how the comparison be- 
tween data and theory develops as the data be- 
come more precise. 

4.3. CP asymmetry in B — > tt + tv~ decay 

The stability of the prediction for the B — > 
■k + tt~ amplitude suggests that the CKM an- 
gle a can be extracted from the time-dependent 
mixing-induced CP asymmetry in this decay 
mode, without using an isospin analysis. Fig- 
ure |^ displays the coefficient S of — sm(ArriB d t) 
in the time-dependent asymmetry as a function of 
sin(2a) for sin(2/3) = 0.75. For some values of S 
there is a two-fold ambiguity (assuming all angles 
are between 0° and 180°). A consistency check of 
the approach could be obtained, in principle, from 
the coefficient of the cos(Ams d i) term, which is 
given by the direct CP asymmetry in this decay. 

5. SUMMARY AND OUTLOOK 

With the recent commissioning of the B fac- 
tories and the planned emphasis on heavy-flavor 
physics in future collider experiments, the role of 
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Figure 4. Mixing-induced CP asymmetry in B — ► 
tt + tt~ decays. The lower band refers to values 
45° < a < 135°, the upper one to a < 45° (right) 
or a > 135° (left). We assume a,/?, 7 G [0, 180°]. 

B decays in providing fundamental tests of the 
Standard Model and potential signatures of New 
Physics will continue to grow. In many cases 
the principal source of systematic uncertainty is 
a theoretical one, namely our inability to quan- 
tify the nonperturbative QCD effects present in 
these decays. This is true, in particular, for al- 
most all measurements of direct CP violation. 
Our work provides a rigorous framework for the 
evaluation of strong-interaction effects for a large 
class of exclusive, two-body nonleptonic decays of 
B mesons. It gives a well-founded field-theoretic 
basis for phenomcnological studies of exclusive 
hadronic B decays and a formal justification for 
the ideas of factorization. 

We believe that the factorization formula (j^) 
and its generalization to decays into two light 
mesons will form a useful basis for future phe- 
nomenological studies of nonleptonic B decays. 
We stress, however, that a considerable amount of 
conceptual work remains to be completed. The- 
oretical investigations along the lines discussed 
here should be pursued with vigor. We are confi- 
dent that, ultimately, this research will result in 
a theory of nonleptonic B decays. 
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